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ABSTRACT In a previous paper we proposed a model to describe a gas of 
pyramidal molecules interacting via dipole-dipole interactions. The inter- 
action modifies the tunneling properties between the classical equilibrium 
configurations of the single molecule and, for sufficiently high pressure, 
the molecules become localized in these classical configurations. The model 
explains quantitatively the shift to zero-frequency of the inversion line ob- 
served upon increase of the pressure in a gas of ammonia or deuterated 
ammonia. Here we analyze further the model especially with respect to 
stability questions. 



1 Introduction 

The behavior of gases of pyramidal molecules, i.e. molecules of the kind 
XY3 like ammonia NH3, has been the object of investigations since the 
early developments of quantum mechanics . In recent times the problem 
has been discussed again in several papers |21|31^|Slini[7] from a stationary 
point of view while in .8. .Q.ilOj a dynamical approach has been attempted. 
For a short historical sketch of the issues involved we refer to [3 E El ■ 

In |H] we have constructed a simplified mean-field model of a gas of 
pyramidal molecules which allows a direct comparison with experimental 
data. Our model predicts, for sufficiently high inter-molecular interactions, 
the presence of two degenerate ground states corresponding to the differ- 
ent localizations of the molecules. This transition to localized states gives 
a reasonable explanation of the experimental results ^2 ^1 . In par- 
ticular, it describes quantitatively, without free parameters, the shift to 
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zero-frequency of the inversion line of NHt, and ND^ on increasing the 
pressure. 

In the present paper we first reconsider our model from the stand point 
of stationary many-body theory clarifying the meaning of the mean-field 
energy levels. We then analyze the mean field states with respect to the 
energetic stability. The conclusions agree with those found in ^U] via a 
dynamical analysis of the same type of model to which dissipation is added. 



2 The model 

We model the gas as a set of two-level quantum systems, that mimic the 
inversion degree of freedom of an isolated molecule, mutually interacting 
via the dipole-dipole electric force. 

The Hamiltonian for the single isolated molecule is assumed of the form 
— -^^tj^, where cr^ is the Pauli matrix with symmetric and antisymmetric 
delocalized tunneling eigenstates and 

Since the rotational degrees of freedom of the single pyramidal molecule 
are faster than the inversion ones, on the time scales of the inversion dy- 
namics set by AE the molecules feel an effective attraction arising from 
the angle averaging of the dipole-dipole interaction at the temperature of 
the experiment Jl]. The localizing effect of the dipole-dipole interaction 
between two molecules i and j can be represented by an interaction term 
of the form —gija^a^ , with gij > 0, where is the Pauli matrix with left 
and right localized eigenstates ipL and (pu 

' (1.2) 
The Hamiltonian for N interacting molecules then reads 
AE ^ 

H = -— ^li«) I2 «)...«) erf (g)...liv 

i=l 
N N 

- ^ ^ li (g) . . . (8) CTf (g) . . . (g) ct| (g) . . . (g) Iat. (1.3) 

2—1 

For a gas of moderate density, we approximate the behavior of the ^ 1 
molecules with the mean-field Hamiltonian 

A E 

h[i,]^-—cj^-G{^,a^^)a\ (1.4) 
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where ?/; is the single-molecule state ((?/', V') = 1) to be determined self- 
consistently by solving the nonlinear eigenvalue problem associated to (|1.4|l . 
The parameter G represents the dipole interaction energy of a single molecule 
with the rest of the gas. This must be identified with a sum over all possible 
molecular distances and all possible dipole orientations calculated with the 
Boltzmann factor at temperature T. Assuming that the equation of state 
for an ideal gas applies, we find 

c^f(f)V.», (1..) 

where Tq — /{AireoerCpkB), £o and Er being the vacuum and relative 
dielectric constants, d the molecular collision diameter and fi the molecular 
electric dipole moment. Note that, at fixed temperature, the mean-field 
interaction constant G increases linearly with the gas pressure P. 



3 Molecular states 

The nonlinear eigenvalue problem associated to H1.4II . namely 

/^[V'mIV'm = M^M (■0^>/.> = 1, (1-6) 

has different solutions depending on the value of the ratio G/AE. If G/ AE < 
i, we have only two solutions corresponding to the delocalized eigenstates 
of an isolated molecule 

= ip+ /ii = -AE/2 (1.7) 
- M2 = +AE/2. (1.8) 

If G/AE > i, there appear also two new solutions 



I AE I AE , ^ 

V'M3 = \/2 + ^ ^+ + V2" 4G A^3 = -G (1.9) 



which in the limit G 3> AE approach the localized states tpL and ipn, 
respectively. Solutions l|1.9|) and (|1.1U|) are termed chiral in the sense that 

The states determined above, are the stationary solutions "ipit) = 
exp(i^t/h)ijjf^ of the time-dependent nonlinear Schrodinger equation 

ih^^m = hi^pim- (1-11) 
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FIGURE 1. Single-molecule energies ei (solid lines) of the four stationary states 
V'Mii * = 1)2,3,4, as a function of the ratio G/AE. The dashed lines are the 
eigenvalues, divided by A*', of the Hamiltonian 11.311 with Qij = G/N and A*' = 12. 



The generic state ip{t) solution of this equation has an associated conserved 
energy given by 

A F C 
f[^]=- — (1.12) 

The value of this functional calculated at the stationary solutions H1.7I1.10|) 
provides the corresponding single-molecule energies Ci — £ [ip^i ] 

ei = -AE/2 

62 = +AE/2 (1.13) 
AE 1 (AE ' 

These energies are plotted in Fig. ^ as a function of the ratio G/AE. The 
state effectively assumed by the molecules in the gas will be that with the 
minimal energy, namely the symmetric delocalized state "0^^ for G/AE < i 
or one of the two degenerate chiral states for G/AE > i. 

The above results imply a bifurcation of the mean-field ground state 
at a critical interaction strength G = AE/2. According to Eq. H1.5|) . this 
transition can be obtained for a given molecular species by increasing the 
gas pressure above the critical value 

9 /T^^ 



Pc. = ^Po ^) , (1.14) 



where Pq = AE/d? . In Tabic [TTTI we report the values of Tq and Pq calcu- 
lated for different pyramidal molecules. 
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AE (cm-i) 


fi (Debye) 


d(A) 


To (Kelvin) 


Pq (atm) 




0.81 


1.47 


4.32 


193.4 


1.97 


ND3 


0.053 


1.47 


4.32 


193.4 


0.13 


PH3 


3.34 X 10-14 


0.57 




29.1 


8.11X 10-1'* 


AsHs 


2.65 X 10-1^ 


0.22 




4.3 


6.44X 10-18 



TABLE 1.1. Measured energy splitting AE, collision diameter d, and electric 
dipole moment fi, for different pyramidal molecules as taken from [151 lltj| . In 
the fourth and fifth columns we report the temperature To and the pressure Po 
evaluated as described in the text. In the case of PH3 and AsH-i the collision 
diameter, not available, is assumed equal to that measured for NH^ and ND^. 
We used Sr ~ 1- 



4 Inversion line 

When a gas of pyramidal molecules which are in the delocalized ground 
state is exposed to an electromagnetic radiation of angular frequency wq ''^ 
AE/h, some molecules can be excited from the state Lp+ to the state (p-. 
For a non-interacting gas this would imply the presence in the absorption 
or emission spectrum of an inversion line of frequency u = AE/h. Due to 
the attractive dipole-dipolc interaction, the value of hD evaluated as the 
energy gap between the many-body first excited level and the ground state 
is decreased with respect to the noninteracting case by an amount of the 
order of G. As shown in Fig. the value of the inversion line frequency 
is actually a function of the number N of molecules and in the limit N ^ 1 
approaches the mean field value 

AE f 2G\^ , , 

'- — ['-ae) ■ ^'-''^ 

According to H1.15|l . the inversion line is obtained only in the range < 
G < AE/2 and its frequency vanishes at G = AE/2. 

In [S] we have compared the mean field theoretical prediction for the in- 
version line with the spectroscopic data available for ammonia [111 I12| and 
deuterated ammonia |13| . In these experiments the absorption coefficient of 
a cell containing NH3 or NDg, gas at room temperature was measured at 
different pressures, i.e. according to H1.5() at different values of the interac- 
tion strength G. The measured frequency v decreases by increasing P and 
vanishes for pressures greater than a critical value. This behavior is very 
well accounted for by the the mean field prediction (|1.15|) . In particular, the 
critical pressure evaluated according to Eq. H1.14|l with no free parameters, 
at T = 300 K is Per = 1-695 atm for NH3 and Per = 0.111 atm for ND3 is 
in very good agreement with the experimental data. 
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FIGURE 2. Inversion line frequency as a function of the ratio G/AE in the mean 
field model (solid line) and obtained from the Hamiltonian 11.31 with gij — G/N 
and iV = 4, 8, and 12. 

5 Energetic stability of the molecular states 

In order to discuss the energetic stability of the mean field molecular states 
found in Section|21we introduce the free energy =£[%()] — niip, ip) ■ The 
stationary solutions of Eq. (|l.ll(l then can be viewed as the critical points 
of the Hamiltonian dynamical system 

Under the effect of a perturbation which dissipates energy, a stationary 
state ip^ will remain stable only if ^[V'/i] is a minimum. Therefore we are 
interested in exploring the nature of the extremal values ^[V'm] of the free 
energy functional. In general, this can be done in terms of the eigenvalues 
and the eigenvectors of the linearization matrix associated to the dynamical 
system (|1.16|l as explained in ^7] for a Gross-Pitaevskii equation. Here, due 
to the simplicity of the model, we can provide a more direct analysis. 

For a variation of the stationary solution -0^ — > -0^ + (50, up to the second 
order in S(j), we have 

^[V^^ + (50] = + (52^[0^, ,50], (1.17) 

where 

A F' 

-G (^(V',,,fT^VM>('50,'T^'^0) + (7/^^,(7^,50) ((50, a>^) 

+ hi;^,a'6^f + i(,50,a>^)2) . (1.18) 
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The variation 6(j) can be taken in the most general form 

Set) = ae''''"(^+ + 6e'''V-, 



(1.19) 



where ip± are the delocalized tunneling eigenstates (|l.l|l and a, 6, 9a, and 
9b arbitrary real parameters with the constraint that + (50, ip^ + S(p) = 
1 + 0((5(/)2). By writing 



(1.20) 



with the real coefficients a^^ and 6^ deduced by Eqs. (|1.7I1.1U|I . the above 
constraint implies 

afj,acos9a + b^bcos9i, ~ 0. (1-21) 

The second variation of the free energy evaluated at the four stationary 
solutions tjjjj,., i = 1,2,3,4, with the condition (|1.21|) gives 



G2b^ (H-cos^e, 
G2a' ( -M-cos^., , , 



(1.22) 
(1.23) 



= G 



2a^ 



AE\ 
^+2g) 



COs2 9a + 2b^ {l~ ^] COS^ (76 



^ I AE , ' 
a\/ 1 - ^ sm T fey 1 + sm 6*6 



(1.24) 



where = 3,4 and the signs =F refer respectively to = 3 and A; = 4. 
We see that, for the state ■0^1 j the variation S'^T is always positive for 
G < and can be negative for G > The variation S^J^ is always 
negative in the case of ipti2- ^'oi' the states ipfj-a and 4' in, which exist only 
for G > the variation S'^f^ is always positive. We conclude that the 
free energy has a single minimum in correspondence of the delocalized state 
■0^j when G < , and two degenerate minima in correspondence of the 
chiral states ■0^3 and -0^^ when G > . The energetic stability analysis is 
summarized in Table lT"^ Note that our results coincide with those reported 
in |1U| where a standard linear stability analysis is performed for the same 
model considered here to which an explicit norm-conserving dissipation is 
added. 



6 Conclusions 



The specific prediction of our model for the critical pressure Per in terms 
of the electric dipole ^ of the molecule, its size d, the splitting AE and 
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-AE 

(52jc-> 
minimum 

saddle point 



+AE 

maximum 
<0 
maximum 



G 
2 



l-(AE 



minimum 



G 
2 



minimum 



TABLE 1.2. Value of the free energy J-[tpti] and sign of its second variation at 
the four extrema ?/>p. , i — 1,2,3, 4. 



the temperature T of the gas, successfully verified in the case of ammonia, 
could be experimentally tested also for other pyramidal gases for which, 
Eqs. 1)1.14(1 and ((1.15|l predict the scaling law 

where 7 = Per X'Y^ I Per XY3 ■ 

Our model applies not only to molecules XY3 but also to their substituted 
derivatives XYWZ . An important difference between the two cases is that 
for XYa the localized states can be obtained one from the other either by 
rotation or by space inversion, while for XYWZ they can be connected 
only by space inversion. This implies that XYWZ molecules at a pressure 
greater than the critical value are chiral and therefore optically active. The 
measurement of the optical activity of pyramidal gases for P > Per would 
allow a direct verification of this prediction. 
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